Functional approach for pairing in finite systems: 

How to define restoration of broken symmetries in Energy 
Density Functional theory ? 

j^>. G. Hupin and D. Lacroix 

Grand Accelerateur National d'lons Lourds (GANIL), CEA/DSM-CNRS/IN2P3, Bvd Henri 
Becquerel, 14076 Caen, France 



o 



M. Bender 



Universite Bordeaux, Centre d'Etudes Nucleaires de Bordeaux Gradignan, UMR5797, F-33175 
Gradignan, France, 

CNRS/IN2P3, Centre d'Etudes Nucleaires de Bordeaux Gradignan, UMR5797, F-33175 
^) , Gradignan, France 



E-mail: hupin@ganil.fr, lacroix@ganil.fr, bender@ganil.fr 



Abstract. The Multi-Reference Energy Density Functional (MR-EDF) approach (also called 
configuration mixing or Generator Coordinate Method), that is commonly used to treat pairing 
£»0 ■ in finite nuclei and project onto particle number, is re-analyzed. It is shown that, under certain 

Q\ ' conditions, the MR-EDF energy can be interpreted as a functional of the one-body density 

00 , matrix of the projected state with good particle number. Based on this observation, we propose a 

l/") . new approach, called Symmetry- Conserving EDF (SC-EDF), where the breaking and restoration 

1^-^ ' of symmetry are accounted for simultaneously. We show, that such an approach is free from 

f***^ , pathologies recently observed in MR-EDF and can be used with a large flexibility on the density 

dependence of the functional. 
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1. Introduction 

In the last decade, configuration mixing has become a standard tool in the Energy Density 
Functional approach both to restore symmetries and to provide a competitive theory for nuclear 
spectroscopy. With the possibility to perform systematic and precise studies, detailed analysis of 
MR-EDF have been made, pointing out some questionable aspects related to the way MR-EDF 
is constructed [U El [31 HJ [5] • Below, we further investigate the MR-EDF applied to particle 
number restoration. The Hamiltonian case is first discussed as well as some specific features 
that will be useful for the EDF case. 

2. Particle number restoration in the Hamiltonian case 

Starting from a two-body Hamiltonian, written here as: 



H = f ij a i a 3 + 4 2^ V ijkl a i a j a l a k , (1) 
ij ijkl 



where v denote the anti-symmetric two-body interaction matrix elements, pairing can be treated 
approximately by minimizing the energy in the space of quasi-particle vacua. For a given quasi- 
particle reference state of the special Bogolyubov type, denoted by \&o), the energy writes 

E S r [$o] = E SR [p,k,k*] = J2 kiPii + VijijPuPjj + \Y, %j K s K iJ » ( 2 ) 

* i,j 1,3 

where p and k denote the normal and anomalous density respectively and where the set of 
single-particle states are those where p is diagonal, i.e. the so-called canonical basis of p. For 
non-zero anomalous density, the U(l) symmetry is explicitly broken and the particle number is 
only imposed in average. 

In a second step, this symmetry can be restored by projecting out the component with N 
particles \^n) = -P^l^o) where P N denotes the particle number projection operator defined 
through [6 

P N = 1_ f 2n dip ^(N-N) (3) 
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The expectation value of any operator O that conserves particle number can then be expressed 

as 



<p N \0\V N ) _ f 2n Jt _ ($o\0\%) 



* I* - / *f (* i* r ^ '^' (4) 
V^nI^n) Jo (<Bo\<B<p) 

where the shorthand 

has been introduced. Here tp denotes the gauge angle, whereas |<E>^) = e lipN |<I>o) refers to the 
state l^o) rotated in gauge space. For O = H, an expression is obtained for the projected 
energy: 

/■27T r -| 

E N [* N ]= dipEsR p 0ip , K 0ip , k v ° Af N (0,tp). (6) 
Jo 1 J 



The energy kernel are defined through 



(®o\H\%) 



Kcd rfV k^ u = x ' 1 Vl (7\ 

E S R[P ,« ,K j - • (7) 

An equivalent form of the projected energy En can directly be found by taking the expectation 
value of H with the projected state, leading to: 

EN[ * N] ~ (*n\*n) ~\ Ui + *U Vijkl ' 



En 



^^l=E%4 + 7E%^<ir (8) 



4 .... 

i j ijkl 



Therefore, we see that the projected energy can be considered either as a functional of the 
transition densities (Eq. @) or as a functional of the one- and two-body density matrices 
(called hereafter OBDM and TBDM) of the projected state (Eq. ©). The two formulations 
are strictly equivalent and can be connected by expressing the density matrices in terms of the 
gauge angles as: 

p% = fj d^M N (0, <p\ = fj *p(p%% - p*>f% + K° k f^)fV N (0, <p) (9) 

Indeed, starting from Eq. (J6j) and combining the pp and the kk terms (Eq. [2]), one can indeed 
directly recognize the components of the TBDM. 



3. Particle number restoration in the EDF case 

The formulation of broken symmetry and its restoration in the nuclear EDF closely follows the 
Hamiltonian case. In a first step, called Single- Reference EDF (SR-EDF), an energy written as 
a functional of an auxiliary quasi-particle state is given. Restricting here to bilinear functional 
based on effective interactions, it can be written as 

£sr[®0] = £SR [p, K,K*]=J2 UiPii + ^ E Vi-ijpiipjj + - E V^jjKijKjj . (10) 

i i,j i,j 

The EDF theory differs from the Hamiltonian case by the appearance of effective vertices v pp 
and v ftK in the mean-field and pairing channels. In particular, these vertices generally (i) differ 
from each other, (ii) might not be fully antisymmetric, (iii) might be density dependent. 

These differences are at the heart of problems occurring when trying to extend SR-EDF theory 
to the Multi-Reference EDF framework. Standard approach to perform MR-EDF consists in 
using equation (jSJ) where Esr is replaced by its functional counterpart Ssr- However, this direct 
mapping has led to inconsistencies: First, MR-EDF can lead to jumps and/or divergences in the 
energy [H O HI [5] . This problem originates in the strategy used to construct MR-EDF and 
requires a specific regularization p2 IH E] that is possible only for specific density dependence 
of the energy. Secondly, the use of density dependent term, that has been so successful in SR- 
EDF is at the origin of technical and conceptual difficulties. Which density should be used in 
the effective kernels still needs to be clarified [7] . Finally, the direct application of Eq. (JBJ) or 
any equivalent for other projection does not insure in the EDF case the proper behavior of the 
energy with respect to the symmetry transformation [8j. None of these problems occur in the 
Hamiltonian case. This is clearly shown by expression p>|h noting that projected densities have 
all the required properties (symmetry transformation, anti-symmetrization, non-divergence and 
regularity...). 

It is therefore quite natural to see if and under which conditions, an equation similar to 
dBJ can be obtained starting from the MR-EDF approach. We have recently shown [9], that 
starting with [a] a bilinear functional, [b] with density independent interaction and [c] a slight 
modification of the renormalization proposed in ref. [3], the MR-EDF obtained by gauge angle 
integration is equivalent to a functional of the projected state densities written as: 

*] = E*^f + \ E i&R&s + \ E ^m R m 

+ lY,^? n ? + ^mn?(l-n»), (11) 

i i 

where denotes the occupation numbers of l^jv). This expression clearly shows that conditions 
[a-c] insure a well behaved functional as far as the particle number projection is concerned. It also 
shows that the MR-EDF can eventually be interpreted as a functional of densities constructed 
from an auxiliary state \fy n) which leads to the functional sequence: 

tfjv => (p n ,R n ) =» £n- (12) 

The direct use of Eq. (jlip instead of Eq. ([6]) is possible and will lead to result identical to MR- 
EDF when the regularization is applicable. This formulation where breaking and restoration 
are accounted for simultaneously will be called hereafter Symmetry-Conserving EDF (SC-EDF). 
SC-EDF can be used with density dependent interaction, including dependence in non-integer 
powers of the density, and can be formulated as a consistent extension of the SR-EDF approach. 
Then, the density entering in the effective vertices becomes the density of the projected state. 
An illustration of the SC-EDF is shown in figure [1] for a functional containing non- integer powers 
of the density. The MR-EDF results where no regularization is possible [8] is also shown as a 
reference. 



Figure 1. (Color online) Particle- 
number restored deformation en- 
ergy surface of 18 O calculated using 
Eq. (jlip when the SLy4 effective 
interaction is used in the particle- 
hole channel. The dashed line 
corresponds to the non-regularized 
MR-EDF result directly obtained 
by gauge angle integration using 
199 points in the discretization (Eq. 

The filled circles correspond 
to the result obtained using the 
Symmetry-Conserved EDF based 
on equation (fTTjh 

4. Summary 

The MR-EDF approach applied to particle number projection has been further analyzed here. 
Guided by the Hamiltonian case, we show that, under certain conditions, this approach can 
be formulated as a functional theory of the projected one- and two-body density matrices. 
This alternative formulation called Symmetry Conserving-EDF can eventually be applied to 
functional where the MR-EDF is ill defined. Note finally that, it has also been shown recently 
that the two-body density matrix of the projected is a functional of the one-body density matrix 
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